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Abstract 

We investigate the leading lattice spacing effects in mesonic two-point correlators computed 
with twisted mass Wilson fermions in the epsilon-regime. By generalizing the procedure already 
introduced for the untwisted Wilson chiral effective theory, we extend the continuum chiral epsilon 
expansion to twisted mass WChPT. Wc define different regimes, depending on the relative power 
counting for the quark masses and the lattice spacing. We explicitly compute, for arbitrary twist 
angle, the leading O(a^) corrections appearing at NLO in the so-called GSM* regime. As in 
untwisted WChPT, we find that in this situation the impact of explicit chiral symmetry breaking 
due to lattice artefacts is strongly suppressed. Of particular interest is the case of maximal twist, 
which corresponds to the setup usually adopted in lattice simulations with twisted mass Wilson 
fermions. The formulae we obtain can be matched to lattice data to extract physical low energy 
couplings, and to estimate systematic uncertainties coming from discretization errors. 



PACS numbers: 11.15.Ha, 12.39.Fe, 12.38.Gc 



I. INTRODUCTION 



A precise matching of results obtained in lattice QCD with the predictions of the chiral 
effective theory is an important test of strong dynamics at low energies. In particular, 
it provides a way to check if chiral symmetry is spontaneously broken according to the 
expected pattern and eventually to extract from first principles the low-energy couplings 
which parametrize the effective theory. 

With many recent calculations with Nj = 2,2 + 1 getting close to the physical point (see 
the plenary talks [H E] presented at the 2009 lattice conference and references therein) this 
goal starts now to be realistic, although a reliable estimation of systematic uncertainties is 
still problematic. In particular, approaching the chiral limit in a finite box implies a detailed 
control over finite-size effects. The chiral effective theory provides information also on the 
volume dependence of physical observables: in the asymptotic region, where M^^L ^ 1, 
volume effects are expected to be exponentially suppressed. For practical purposes, for 
instance for the extraction of the pion decay constant, the empirical criterion Mt^L > 4 
seems to be necessary in order to keep effects well below the statistical error typically 
present. 

An alternative approach is to study QCD in a different kinematic corner, namely the e- 
regime |31ll]: here finite- volume effects are polynomial instead of exponentially suppressed, 
and one can exploit the finite-size scaling properties of given observables in order to extract 
information about infinite-volume quantities. Since the chiral expansion obeys a different 
power counting with respect to the infinite volume, higher order corrections will be different: 
the matching of the chiral effective theory with lattice QCD will then provide low energy 
couplings which will be affected by different systematic uncertainties. A general agreement 
among those independent determinations is a good check on the validity of the approach. 

Many quenched simulations in the e-regime have been performed [514T3]. using Dirac 
operators which satisfy the Ginsparg-Wilson relation [l^. More recently, e-regime calcula- 
tions with Ginsparg-Wilson fermions have been carried out also in the dynamical case, with 
Nf = 2 [T3H2n] and Nf = 2 + 1 [211 [22] • Since the Ginsparg-Wilson relation ensures exact 
chiral symmetry at finite lattice spacing |23] , it guarantees many theoretical advantages, for 
instance the possibility to reach arbitrarily small quark masses and a continuum-like renor- 
malization pattern. The price to pay is the high computational cost, which makes these 
simulations very challenging. In particular, approaching the continuum limit or exploring a 
broad range of physical volumes requires very big efforts. 

Recently it has been realized that simulations in the e-regime are feasible also for Wilson- 
type fermions. In pH - [26] first results obtained with Wilson twisted mass fermions have 
been presented. The use of a PHMC algorithm combined with an exact reweighting of 
a few low modes of the lattice operator turned out to be an important ingredient in this 
study. Analogously, in [23 EH] a reweighting algorithm has been proposed and successfully 
applied to simulate nHYP improved Wilson fermions in the e-regime (see also ref. [29]). 
In both cases continuum chiral perturbation theory (ChPT) describes the lattice data very 
well, although chiral symmetry is explicitly broken for Wilson-type fermions. Even though a 
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scaling study would be necessary to systematically investigate lattice artifacts, this suggests 
that the impact of chiral symmetry breaking is mild and that it can be legitimate to match 
lattice results with the expressions of continuum ChPT. 

We address this issue by means of Wilson Chiral perturbation theory (WChPT) [301 EI] , 
which can be generalized to the twisted mass case [32] - [37j . In Refs. [Ml [39] we extended 
untwisted WChPT to the e-regime.^ 

The relevant issue is the relative power counting of the quark mass and the lattice spacing 
(in units of Aqcd)- It turns out that for m ~ qAqqj-, (GSM regime) the explicit breaking 
of chiral symmetry is still dominated by the quark mass, and lattice artifacts are highly 
suppressed. For mesonic two-point functions, the lattice spacing corrections start to appear 
at NNLO. On the other hand, if m ~ a^AQ,^^, (Aoki regime), lattice artifacts compete with 
the quark mass, and corrections are substantial since they contribute already at LO. 

In Refs. [38], [39] we also introduced an intermediate regime (GSM*), where discretization 
effects appear at NLO, and for this case we computed the leading corrections for several 
correlators. An important observation is that in this intermediate regime only one additional 
low energy coupling appears, namely C2 [30J. In this paper we extend this study to the 
twisted mass case. While there will be many features in common to the untwisted case, 
some new aspects arise and will be discussed. 



The paper is organized as follows: in section [IT] we define the chiral Lagrangian of twisted 
WChPT, currents and densities, and we recall the main related properties and definitions in 
an infinite volume; in section III we define the effective theory in the e-regime, we discuss the 
power counting and the role of the vacuum; in section IV we compute leading corrections in 
the GSM* regime for several 2-point correlators and we give numerical estimates for these 
corrections. Finally, we draw our conclusions in section |V| 



II. CHIRAL PERTURBATION THEORY FOR TWISTED MASS WILSON 
FERMIONS 

A. The chiral Lagrangian 

Correlation functions computed with lattice simulations are affected by discretization 
errors, which can be analyzed using effective field theory. To obtain the correct form of the 
effective Lagrangian one proceeds in two steps [3^. First, one matches the lattice action used 
in the simulations with the appropriate Symanzik effective action. The Symanzik action is 
subsequently matched to a chiral Lagrangian which contains the standard continuum terms 
and appropriate additional operators that transform under chiral symmetry as the operators 
of the Symanzik effective theory. These additional operators describe the effects of the 



^ Recently, WChPT in the e-regime has been adopted also to study the spectral density of the Wilson Dirac 
Operator at fixed topology [iO] . 
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nonzero lattice spacing a? 

In this paper we are interested in lattice actions with Wilson twisted mass fermions |13] . 
These actions are sometimes called Wilson-type fermions because they represent a simple 
generalization of the standard Wilson action. Their explicit form will not be needed in 
this paper (see ref. [Hj for a review). We just recall that in our analysis we consider the 
lattice action with A^f = 2 degenerate flavours, and bare mass parameters mg and /iq. The 
untwisted quark mass = Zm{niQ — nicr) and the twisted quark mass /ir = Zp^fi^ are 
renormalized with renormalization factors computed in a mass independent scheme, and rricr 
denotes the critical mass. Wilson twisted mass lattice QCD in the continuum limit is then 
equivalent to QCD [13] with a physical quark mass mp = ^/m'^ + /x^. From now on we will 
drop the subscript R and all the quark masses are considered, unless specified differently, as 
renormalized in a mass independent scheme. 

The chiral Lagrangian in twisted mass WChPT is essentially the same as for the untwisted 
case [33 ll35ti371 H5] . The only difference is a mass matrix that contains, besides an untwisted 
mass m, the twisted mass /i. 

In WChPT there are two sources of explicit chiral symmetry breaking, the quark masses 
m, fi and the lattice spacing a. The power counting is determined by the relative size of these 
parameters. The literature [Ml distinguishes two regimes with different power counting: 

(i) the generically small quark mass (GSM) regime where the quark mass is ~ aA'^^^i 

(ii) the Aoki regime where the quark mass is ~ o^Aqq^)- Depending on the particular 
regime, the LO Lagrangian differs. Explicitly, in the GSM regime at leading order the chiral 
Lagrangian reads 

£2 = ^Trid^Ud^U^) - ^TiiM^U + U^M) . (1) 

F and S are the familiar low energy couplings (LECs) appearing at LO in the continuum 
chiral Lagrangian and U contains the pion fields in the usual way (see below). The mass 
matrix Ai is defined by (a" denotes the Pauli matrices) 

M = (ml + ifia^). (2) 

The field U and the matrix Ai are written in the so called twisted basis. There is no 0(a) 
contribution in the chiral Lagrangian C2, since we have absorbed it in the definition of m, 
which therefore represents the so-called shifted mass [30].^ 

Eq. (|2| can be used to define the polar mass mp and a twist angle coq by 

M = mpe^'^o"' , (3) 

or, directly in terms of m and fi: 

mp = A/m^ + , tan Wo = — . (4) 

m 



^ For introductory lecture notes see Ref. [HI 02] . 

^ The shifted mass is sometimes denoted by m'. Here we drop the prime since we exclusively work with the 
shifted mass and we will never need the original mass parameter. 
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The subscript '0' serves as a reminder that the angle uq relates the two mass parameters in 
the chiral Lagrangian.^ We emphasize that uq is related to but not identical with the twist 



angle u that we define in section II C 



We can easily go the so-called physical basis, where the quark mass matrix takes the 
standard form proportional to the identity, Ai = mpl, by performing the following non- 
anomalous field transformation: 

U = WUW with W = exp (^^^^^ • (5) 

Here U denotes the field in the physical basis. 

If we move to the Aoki regime we have to add to the Lagrangian Q the O (a ) corrections 

[161 HZ] 



^2 



<5£,2 = — C2a2(Tr(f/ + f/t))% (6) 

since these contribute already at LO in this regime. Therefore, we expect the corrections due 
to a nonzero lattice spacing to be much more pronounced in the Aoki regime. Additional 
chiral logarithms proportional to appear in one-loop results of various observables like 
the pion mass and pion scattering lengths [ITJIIH]- Moreover, in infinite volume non-trivial 
phase transitions become relevant, and the sign of the LEG C2 plays a decisive role for the 
phase diagram of the theory [30] . 

The NLO terms linear in a can be found in Refs. ^ST[ SB] , but these terms will not be 
needed in the following. 



B. Currents and densities 

The currents and densities in the effective theory can be computed by adding appropriate 
source terms to the partition function of the theory and differentiating the resulting generat- 
ing functional with respect to the sources [36lll9]. Alternatively, the currents and densities 
can be obtained using a standard spurion analysis, as described in Both methods have 
been discussed in detail in the literature. Here we simply summarize what is relevant for 
our computation with Wilson twisted mass fermions. 

We assume the following definitions for the currents and densities at the quark level: 

S°{x) = ^{x)^{x), P\x) = z^(x)75TXx), (7) 

A-^ix) = i^{xh,j,T^^P{x), V^{x)=ti^{xh,T^^P{x). (8) 

T", a = 1, . . . , Nj — 1 are the Hermitean SU{Nf) generators satisfying the property 
Tr(T'^T'') = 6°'^/2. For Nf = 2, = a°'/2. The corresponding currents and densities 



^ Since we use the shifted mass to parametrize the chiral Lagrangian, the angle luq defined in Q does not 
coincide with the one in Ref. [55] . 
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at LO in the chiral effective theory read 

= -^Tr [U{x) + U^{x)] , (9) 

= i^Tr [T^iUix) - U\x))] , (10) 

^^(a;) = * y'^' [r'^(f/(x)t9^f/(x) - t/(x)9^f/t(x))] , (11) 

V;(x) = -^^Tr [T"(f/(x)t9,,f/(x) + f/(x)9^f/t(x))] . (12) 

These are the famihar expressions from LO continuum ChPT. The leading corrections of 
0(a) can be found in [SSIHU] . They are of higher order and we will not need them explicitly 
in the following. 



C. The PCAC mass and the twist angle 

There are various ways to define a twist angle. A popular way that has often been used 
is by the ratio 

tancj = — - — . (13) 

"^PCAC 

The denominator is given by the PCAC mass, which is an observable, instead of the shifted 
mass m. The former is defined as usual by 

(a^A°(x)P"(y)) 

^^^^^ - 2{P^{x)P^{y)) ■ ^^^^ 

In the classical continuum limit the PCAC mass is equal to m, and therefore u = uo- For 
nonzero lattice spacings this is true only at LO, but violated at higher order, as we will 



derive in section [TV B[ Hence, usually we have oj ^ uq. 

Quite generally, mpcAc is a function of m, /i, the lattice spacing a and the volume V. In 
principle this function can be "inverted" to obtain the shifted mass m as function of the 
other parameters, 

"^PCAC = '^PCAc('^,/^,a, V") =^ m = m{mpcAc,l^,a,V) . (15) 

In practice this relation can be computed perturbatively in the chiral expansion. Once it 
is known it allows us to express other observables as functions of the measurable quantity 
'^PCAC instead of m. This will be the way we present our results in the latter sections. 
A particularly interesting value for the twist angle is maximal twist, u = 7r/2, where 



automatic 0(a) improvement is guaranteed [50]. According to our definition (13) maximal 
twist is given by a vanishing PCAC mass. 
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III. EFFECTIVE THEORY IN THE EPSILON REGIME 



A. Epsilon expansion and power counting 

The discussion on the power counting in the previous section assumed an infinite volume. 
Finite-size effects due to a finite volume V = TL^ with L,T ^ 1/ Aqcd can be systematically 
studied within ChPT [3l HI |5T] . If the pion Compton wavelength is much smaller than the 
size of the volume, M^^L ^ 1, finite- volume effects can be treated in the chiral effective 
theory by adopting the standard p-expansion, where the inverse box extensions are treated 
as small expansion parameters of the same order as the typical momenta: 1/L, 1/T ~ 0{p). 
For asymptotically large volumes the finite- volume corrections are exponentially suppressed 
by factors exp(— M^L). 

On the other hand, for pion masses such that M^L < 1 one probes the so-called e- 
regime of QCD [52H51]. In this regime the pion zero mode contribution to the 

pion propagator will eventually diverge in the chiral limit. Hence its contribution cannot be 
treated perturbatively but has to be computed exactly. This is achieved by a reordering of 
the chiral expansion by means of summing up all Feynman graphs with an arbitrary number 
of zero modes propagators. In the e-expansion one parametrizes the chiral field U according 
to 

U{x) = Uoexp(^^C{xK^ , (16) 

where the constant Uq G SU (Nf) represents the collective zero-mode. The nonzero modes 
parametrized by are still treated perturbatively. These satisfy the condition 

/ rf^xC(x) = 0, (17) 
Jv 

since the constant mode has been separated. The e-expansion is now defined by using the 
counting rules 

M2~0(e^), 1/T, 9^ ~ 0(e), T ~ 0(e). (18) 

With these counting rules the product M^V counts as e°, just as its inverse. Consequently, 
all Feynman graphs that exclusively involve zero-mode propagators count as 0(1) and are 
unsuppressed. 

Once the counting for the pion mass is fixed, it determines the counting of the quark mass. 
In continuum ChPT the tree- level result = 2Bm, with B = H/F"^, fixes the counting 
m ~ 0(e'^). The same line of argument has been applied to WChPT for untwisted masses 
m, and one obtains the same counting [38l |39]. Here we use it to deduce the counting for 
the polar mass mp and the twisted mass fi. The continuum tree-level pion mass for twisted 
mass QCD is = 2Bm-p. Counting the pion mass squared as O(e^) we immediately find 

mp ~0(e^). (19) 

This is the result that we would naively expect. With Wilson twisted mass fermions mp 
plays the role of the physical quark mass, so the counting should be as for m in the untwisted 
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case.^ In terms of m and fi this corresponds to have either both masses of O(e^) or at least 
one of them of O(e^) and the other of even higher order (in case one of the two masses is 
significantly smaller than the other). 

In our computations we will keep both masses to be of O(e^). This choice allows us to 
keep the computation general and our final formula, expressed in terms of rescaled masses 
Zm = m'EV and = /iSV^, are valid for arbitrary twist angles. In particular, our final 
results will account for the special cases u = tt/2 (maximal twist) as well as for a; = 0, 
where we reproduce the results for standard Wilson fermions obtained in [381 EH] • 

The counting rule for the lattice spacing a is now easily fixed. Quite generally, as in 
infinite volume, the power counting is determined by the relative size of a and the quark 
mass. Since the counting of mp is fixed we obtain the counting for a. The arguments are 
just as in WChPT with untwisted masses [SSIES], in particular, we carry over the definitions 
for three different regimes: 



Depending on the relative size of a and mp one of these counting rules is applicable. The 
GSM and the Aoki regimes have also been introduced in infinite volume WChPT. The e- 
expansion allows the introduction of the intermediate GSM* regime [39], which defines the 
"transition region" between the other two regimes |38j . 

In the e-regime the topological charge plays a relevant role [53| and predictions in sectors 
of fixed topology can be given in the chiral effective theory. Since for Wilson-like fermions 
there is no unambiguous definition of the topological charge at non-zero lattice spacing, we 
compute correlators where all sectors have been summed up. 

B. Vacuum state and epsilon regime 

A slightly unusual feature in twisted mass WChPT is the non-trivial ground state Uv, 
that is determined by a gap equation [35l |37|. Uy depends on the parameters m,/i and a, 
and this dependence affects observables when computed perturbatively. The reason is that 
Uv enters the calculation if we compute correlation functions perturbatively by a saddle 
point expansion of the path integral around Uy- In the following we want to argue that Uy 
is no longer needed in e-regime calculations where one integrates exactly over the collective 
constant mode. 

In chiral perturbation theory (continuum or on the lattice, with a twisted or untwisted 



^ Notice that in the Aoki regime the situation can be more comphcated, since the standard LO relation 
between the pion mass and the polar mass is modified by O(a^) corrections. 



GSM regime : 
GSM* regime : 
Aoki regime : 



a ~ O(e^) , 




(20) 
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mass) we are interested in correlation functions defined by a functional integral, 

{O) = ^ J V[U] 0[f/]e-^-[^l , (21) 

where 5"^^ is the effective action and O an effective (local) operator at a given order. VlU] 
denotes the measure for the path integral that needs to be properly defined [56l [57] . 

Suppose the effective action assumes its minimum for the constant field configuration 



U = Uv and the integrand in (21 ) is strongly peaked around it. In this case we perform the 



standard saddle point expansion around Uv- We expand the field U by the familiar ansatz 

U{x) = Uvexp{i7Ta{x)ajF), (22) 

and the measure is given by the formal product measure 

V[U] = V[n], (23) 

leading to standard Gaussian integrals involving the propagator for the pion fields. 

Alternatively, we may parametrize the U field by isolating the collective zero-mode field 
as it is done in e-regime calculations. In this case we write 

U{x) = UvUoexpiiUxW/F), (24) 

where the integration over Uo is done exactly. The measure for this parametrization then 
reads [H |52] 

V[U] = d[f/o]P[e][l-5(0]- (25) 

The measure factorizes and the correction [1 — i?(,^)] can be exponentiated to give an effective 
action ^measf^]- In Appendix [A| we give the exphcit expression of B{^) at O(e^). The main 
point here is that the measure d[t/o] for the constant mode is the standard Haar measure 
over the group manifold SU{Nf). Since it is left-invariant we can use d[Uo] = d[UvUo] in 
computing the path integral. This implies, since we integrate over all constant fields, that 



the particular field Uy is irrelevant, and instead of (24) we can directly parametrize the 
fields according to 

Uix) = Uoexp{tUx)aa/F), (26) 



just as one does in calculations with an untwisted mass term [3S1 ES]- To summarize: In 
contrast to twisted mass WChPT in infinite volume (or in the p-regime) the ground state 
and the gap equation do not play a special role in the epsilon regime. 

C. Epsilon expansion of correlation functions 

In the following sections we mainly use the notation of Ref. |3pJ. In appendixjCjwe briefly 
discuss how the final formulae should be written using the notation of Ref. [38] . 
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The discussion of the e-expansion of correlation functions can be carried over from Ref. 
[38| [39] . It is based entirely on dimensional arguments and once m is replaced by mp the 
entire discussion holds true for the twisted mass case. For this reason we do not repeat the 
arguments here but simply summarize the main results. 

Correlation functions {Oi(x)02{y)) = (O1O2) (for notational simplicity we suppress the 
dependence on x, y) in WChPT are written as the sum of the corresponding continuum 
correlator plus a correction stemming from the nonzero lattice spacing, 

WChPT — (Oi02)ct + 5(Oi02). (27) 

The correction 6{Oi02) receives contributions from both the effective action and the effective 
operator. 

This correction is proportional to powers of the lattice spacing. At which order it con- 



tributes depends on the regime, cf. eq. (20). 

In the GSM regime the correction starts with higher than the continuum contribution. 
In other words, the lattice spacing first affects the correlators at NNLO. Working to NLO 
one can ignore the correction and the continuum results are the appropriate ones. 

In the GSM* regime the lattice spacing corrections enter at NLO. However, at this order 
only the O(a^) correction proportional to C2 contributes. The corrections linear in a, stem- 
ming from the corrections in both the action and the effective operators, are suppressed by 
one more power of e.^ 

Most pronounced are the corrections in the Aoki regime, where they contribute already 
to LO. In addition, the O(a^) correction in the chiral Lagrangian cannot be completely 
expanded, it provides a zero-mode contribution of order e° that has to be treated exactly. 
As a result, the integrals over the constant mode are no longer the standard Bessel functions 
that one usually encounters in e-regime calculations. 

Notwithstanding the complications in the Aoki regime, the main conclusion one can draw 
is that the lattice spacing corrections are typically suppressed, in the GSM regime to NNLO. 
This suppression of the lattice spacing corrections is one of the main reasons for the belief 
that Wilson fermions (twisted or not) are still a good choice for e-regime simulations despite 
their explicit chiral symmetry breaking. 



IV. LEADING CORRECTION IN THE GSM* REGIME 
A. Basic definitions 

Lattice spacing corrections to correlation functions enter at NLO in the GSM* regime. As 
mentioned in the previous section, only the O(a^) term in the effective action contributes. 



® These corrections of 0(am) for Wilson fermions have been computed in Ref. [35] for the pseudoscalar and 
scalar two point functions. 
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and the NLO correction explicitly reads 



6{0,ix)0,iy)) = -{0'^^,ix)0;^^,,iy)5Sa2) + {0'^^,ix)0;^^,,iy)){6Sa:) . (28) 

GSM*, NLO 



The superscript "LO" in (28) refers to leading order in the e-expansion. In the twisted basis 



the O(a^) term takes the same form as for standard Wilson fermions 

SS^. = {-mUo + U^o)Y, (29) 
where we introduced the dimensionless quantity 

p = F^C2a'^V. (30) 



The angled brackets in (28) stand for the functional integral over the non-constant fields 
^"(x) and the constant mode Uq. The integrals over the first ones are done perturbatively. 
This part is completely analogous to the untwisted case in Ref. El], and we refer to 
appendix |A] for a collection of useful properties of the propagators. 

The integral over the constant mode has to be done exactly, and here differences appear 
because of the twisted mass term /i. The integrals we encounter in the twisted basis are of 
the type 



(giUo)) = ^ I d[Uo] g{Uo) e^^[^o+^oV^^T.[.3(^„_^r)] ^ ^^^^ 



^0 JSU{2) 

where Zq is the partition function, obtained with g{Uo) = 1. The parameters and are 
defined as'' 

Zm = mYy , z^ = fiTiV . (32) 



For = the integral (31 ) reduces to the standard one for an untwisted mass term. In this 
case it leads to expressions involving modified Bessel functions In{x) with integer index n. 

The same is true for z^ ^ 0, although the integral looks superficially rather different. This 
is immediately seen after performing a field redefinition to the physical basis. Performing 
the axial rotation ([s]) and using the invariance of the Haar measure, d[Uo] = d[WUoW], eq. 



(31) can be written as 

{giUo)) = ^ I d[f7o] ^7(W-f^W^) ei^=^[^°+^'l ^ {9{WU,W))^^,, , (33) 

^0 JSU{2) 

where z in the exponent is now given by 



mpSl^ = Jzl + zl- (34) 



^ Zra is the standard combination familiar from continuum ChPT, where it is usuahy denoted by /x. In 
order to avoid confusion with the twisted mass we had to change this notation. 
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The Boltzmann factor in the integrand assumes now the standard form with the polar mass 



mp entering the exponent, hence the index 'phys' on the r.h.s of eq. (33 ). The representation 



(33) is particularly useful for doing actual calculations, since many results for integrals with 
untwisted masses can be taken over to the twisted mass case. 

Note that twisted mass Wilson fermions break the SU(2) isospin symmetry to a residual 
U(l). This implies that completeness relations for SU(A^/) generators, which are repeatedly 
used in standard e- regime calculations [52], cannot be apphed in our case. This poses a 
slight computational nuisance but no serious difficulty. In Appendix |D] we give an example 
for the computation of zero mode integrals in presence of isospin breaking. 



We calculated the correction (28) for a variety of mesonic correlation functions. For 
the presentation of our results we find it useful to introduce the following notation. First, 
translation invariance allows us to write 



{X\x)Y\y)) = Cfy{x-y), (35) 

where X" and F° represent one of the densities or currents listed in ([9]) - (12). We suppress 
the spacetime index in the currents and only make the isospin index explicit. The vector 



and axial correlators we consider in section IV C 2 refer to the temporal component of the 
currents. 

In the GSM* regime, CxY{^~y) can be written through NLO as the sum of the continuum 
correlator and a correction proportional to a^, 

CI'y ix-y) = C^xY,c, {x-y) + C^xY,a^ (x-y). (36) 

The continuum correlator (for generic Nf) at NLO can be found in the literature [52] (see also 
Refs. [HH] and [SH])- Since we present our final results in the twisted basis, for the reader's 
convenience we collect a few relevant formulae and the expressions for the continuum parts 
of the correlation functions in the twisted basis in appendix [B| 

For the matching with numerical results obtained in lattice simulations one is often in- 
terested in the correlation function integrated over the spatial components, 

C^yit) = [ d'xC%{x -y) = Cty,At) + C^y^a^ , (37) 

To the order we are working here the correction C'^ya'^ is independent of t and only shifts 
the constant part of the continuum result. 



B. The PCAC mass in the GSM* regime 



The first observable we compute is the PCAC mass defined in (14). This allows us to 



express the results for other correlators as a function of mpcAc instead of m (or, equivalently, 
as a function of u instead oi uq). 



In the chiral effective theory we write the numerator in (14) as (no summation over the 
fiavor index a, with a restricted to 1,2) 

{d,A^^{x)V''m = CaAPix) , 

CdAp{x) = C9AP,ct{x) + CoAP,a2ix) , (38) 
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and similarly for {V°'{x)V°'{0)). To leading order in the e-expansion we find 

S l2{2z) 



CdAP,ct (x) 

Cpp^x) 



V h{2z 

h{2z) 
2 zh{2z) 



cos Cdo : 



where z is defined in (34). Hence, for the PCAC mass we find the expected result 

m-p cos (jJq 



m 



LO 

PCAC 



m , 



(39) 
(40) 

(41) 



Defining maximal twist by a vanishing PCAC mass is therefore equivalent to m = 0, at least 
to this order in the chiral expansion. 



The leading correction in the GSM* regime is given by (28). The NLO result for the 
PCAC mass is then found to be 



NLO 
PCAC 



m-p cos Wo (1 + pAm) = m (1 + pA^ 
2 h{2z) 



zl2{2z) 



(42) 
(43) 



The main observation we can make is that the PCAC mass is still proportional to m. 
Consequently, maximal twist, given by a vanishing PCAC mass, is still equivalent to m = 0. 



This is perhaps better seen if we reformulate (43) in terms of the twist angles 



tana; = (1 — pAm) tanwo . 



(44) 



This implies that w = 7r/2 is equivalent to = 7r/2 at NLO. This is another way of saying 
that ?T^pcAc vanishes if m = 0. Note that u = Uq only for angles ±7r/2 and 0. For all other 



values full use of (44) has to be made. 



Eq. (42) can be inverted to obtain m as a function of mpcAC, 

"^PCAC [1 - pAm] , 



m 



(45) 



where in the correction A^ we can use the LO result (41) and replace m by mpcAc- This 



is the result we already anticipated in eq. (15). In the following we will use it to express 



correlators as functions of mpcAc instead of m. 



C. Results 

Our results are presented in terms of 

Zm = mpcAcSV, 



pSV. 



(46) 



For Zm we keep the same symbol as in eq. (|32|), but we now substitute m with mpcAc given 

(47) 



in eq. ( 45 ) . In order to make our formulae more readable we also use 



^fi ~^ ^ni- 



ls 



1. Scalar and pseudoscalar correlators 



It is convenient to write the O(a^) correction to the time correlator in eq. (37) in the 
form 



^PP,a'2 



ab 
PP^ 



c 



00 

SS,a^ 



00 
55- 



For the pseudoscalar correlator we obtain 

1 



11,22 



2zHi{2zYh{2z 
- z + llzj) h (2^) h {^zf + 

+ 2 {-2Zra^ + 'izj + (1 - 2zJ)) h {^zf h {^z) + 

+ z + 5zJ) h (2^)' 



(48) 



(49) 



A 33 
LApp 



4 {-zj + Az^^zJ + bz^^zj) h {2zf + 



2z7/^(2z)2/2(2z 
+ z (5^/ + 42z„2^/ - llzj) h {2z) h {2zf + 

+ 2 (-2z„6 + 2zJ + (2^„2 - 1) + (8^„^ - 2^„2)) h {^zf h (2z) + 



(50) 



+ z (-3z/ - Uzjz^^ + 5zJ) h (2x 
The Wilson untwisted case corresponds to = 0. In this limit we obtain 

^ab 



Aab 
L\pp 



-0 2zlh{2z^Yl2{2zm) 



Azih{2z^Y - llzmh{2zjh{2zmy + (51) 



+ 2 (3 - 2zl) h{2z^fh{2z^) + bz„M2zmfy 

which reproduces the result reported in [32] (eq. 4.46) and [M] (eq. 4.4). 

The most interesting case is maximal twist, which is obtained for Zm = 0. In this case 
the results simplify to 



11,22 



A 33 
^PP 



-zM2z^f + 2h{2z^)h{2zyj + ^^/2(2^^)^ 

2z//i(2z^)2 
5z^/i(2z^)^ - 10/2(2^^)/i(2z^) - 3^^/2(2^^)^ 
2^^/i(2;^,)2 



(52) 
(53) 



11 22 

Notice that both corrections App and App are finite in the limit z^ — )• 0. 



The O(a^) corrections for the scalar singlet and pseudoscalar correlators are related by 



AGO 

^ss 



a=l 



aa ri 

pp — U, 



(54) 
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which is vahd for generic twist angle, at least to this order in the chiral expansion. This 
result leads to a generalization of what has already been found for untwisted masses [38] : 
The combination 



+ J2CMt) (55) 



a=l 

of correlation functions is free from Ofa^) corrections. 



2. Axial and vector correlators 



For the O(a^) correction in eq. (37) to the time-component axial and vector correlators 
we find 

where A"^^ is given by 



xY,a^ - -P^AxV, X,Y = A,V. 



(56) 



A 



11,22 
AA 



-^—^\ 4 {-z^' + z^^zj + 2z^''zJ) h {2zf + 



(57) 



+ z (-2z/ - 15zjz^^ + llzj) h (2^) h {2zf + 

- 2 [-2zJ + ^Zm^ + z^^ [AzJ - 2) + z^^ {2zJ + zj)) h {2zf h {2z) + 

+ Z [z^'^ + AZm^Z^^ - 5Zm^) h (2z)' 



A 33 
^AA 



Azjz^h i2zf + z {z^^ + llzj) h {2z) h {2zf + (58) 



z^h{2zYl2{2z) 

+ 2 {{2zJ - 1) z^^ + zj {2zJ - 3)) J2 {2zf h {2z) - z {z^^ + 5zJ) h {2z 



Also in this case one can verify that for = one obtains the Wilson untwisted formula 



^AA 



w=0 



-2Afp 



ui=0 



At maximal twist we obtain 
A 11,22 _ 

^AA,uj=tt/2 ~ 



(59) 

(60) 
(61) 

The O(a^) correction for the vector correlator is, up to a sign, the same as for the axial 
vector correlator. 



-2z^h{2z^f + Al2{2z^)h{2z^) + ;^^/2(2^^)^ 
^lhi2z,r 



,33 _ z^h{2z^f- 2h{2z^)h{2z^) - z^h{2z^f 

^AA,uj=-k/2 — 



A a6 \o,b 



AA,a'^ 



(62) 
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This identity holds for generic twist angles and generalizes the result for untwisted masses 
[39] . An obvious consequence is that for correlation functions of right- and left-handed 
currents ^ ji= |(V^±A^) the leading O(a^) corrections cancel. This follows from the fact 
that these currents do not contain zero modes at LO, hence the connected and disconnected 



contributions in eq. (28) cancel among each others. 



D. Numerical estimates 



For the correlators considered, the leading O(a^) correction in the GSM* regime is just 
a shift of the constant part. In order to get estimates for the size of these corrections in 
typical present-day simulations we look at the ratios 



jDab 
rLxY 



(63) 



These ratios are the relative shift of the correlators at the midpoint t = T/2. 

Approximate values for the parameters entering the correlators are taken from the sim- 
ulations of the ETM collaboration [261 [59]. We use F = 90 MeV and a = 0.063 fm. For 
simplicity we assume a hypercubic lattice with Nt = Nl = 24, which corresponds to a 
box size L = 1.512 fm. This is rather small and the e-expansion might not converge well, 
but here we are interested only in an estimate about the order of magnitude for the lattice 
spacing corrections. 

The coefficient C2 is estimated from the pion mass splitting together with the LO ChPT 
prediction — 2c2a^ = — m^o |36]. The data for the charged and neutral pion masses 
in Ref. [59] translates into \c2\ ~ (600 MeV)^. The error, however, is quite large because 
of the large statistical uncertainty in the determination of the neutral pion mass.^ |c2| ~ 
(600 MeV)'^ implies p ~ 0.75. Although this is slightly large it is smaller than 1 and we may 
still count this as O(e^), as we should in the GSM* regime. 

Figure [I] shows -Rpp^, -Rpp, Ryy'^ and R'yy for maximal twist {zm = 0) and z = values 
in the e- regime. For = 1.0 we find values less than 2.5 percent, decreasing to less than 
2 percent for = 2.5. The ratios for fiavor index a = 6 = 1, 2 and a = b = 3 assume the 
same value for vanishing z^, since for vanishing twisted mass we restore isospin symmetry. 
At z^ = the corrections are maximal but for small values 2;^ we eventually enter the Aoki 
regime and our formulae cease to be valid. 

The curves in fig. [T] look qualitatively very similar to the ones for untwisted masses 
shown in Ref. [39], although the decrease of the ratios for growing z^ is slightly faster in the 
untwisted case. 

Figure [2] shows the ratios involving the axial vector current, R]^^'^ and R^j^- The ratio 
R^^ shows a behavior similar to the ratios plotted in figure [I] In contrast, R]lf^ increases 

® Note that the vahie for C2 is not universal but depends on ah the details of the lattice action chosen in 
the simulation. An analysis [^01 of quenched twisted mass lattice data led to the value C2 ~ (300 MeV)^. 
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FIG. 1: Ratios for the pseudoscalar and vector current correlators as a function of z 



with increasing z^, up to about 6 percent for 2;^ = 2.5. The reason for this somewhat odd 
feature is not that the O(a^) corrections are larger for this particular correlator. The origin 
for the increase in -R^^jf ^ is the continuum correlator in the denominator of the ratio. C^a^aIx 
deceases much more rapidly with increasing than ^2 in the numerator, leading to an 
increasing ratio -R^^jf^- In fact, the result for B^^'^ will eventually diverge for large 2;^ where 
the continuum correlator has a zero. Notice that this happens in a region where ^ ^ 1 (at 
fixed volume), which is not expected to be in the domain of validity of the e-expansion. Fig. 
[3] shows directly the correlators Cj^^^^(T/2) and C^^'^l^iT jl), both divided by {—F'^/T) in 
order to get dimensionless quantities. Obviously, the correlators are well behaved and 0{a?) 
correction gives a small and almost constant shift of the continuum result. 

We conclude that, for our choice of parameters, the O(a^) corrections to the correlators 
are at the few percent level, a small and probably negligible correction. 



V. CONCLUDING REMARKS 



We have extended the framework of the Nf = 2 Wilson chiral effective theory in the e- 
regime to the case of a twisted mass. By keeping the same power counting for the untwisted 
mass m and for the twisted mass /i, we have defined different regimes along the same lines 
as in the pure Wilson case [381 El]. For quark masses of the order aA^^-^j^ (GSM regime), the 
explicit breaking of chiral symmetry induced by the lattice spacing is strongly suppressed, 
and lattice artefacts appear only at NNLO in the epsilon expansion. On the other hand, if 
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FIG. 2: Ratios for the axial vector correlators as a function of z 



the quark masses are of order a^Aq^j-, (Aoki regime), discretization effects appear already 
at LO. In this paper we have focused on the intermediate (GSM*) regime, where lattice 
artefacts start to contribute at NLO, which is the order at which the matching between 
lattice data and the chiral effective theory is usually performed. We have computed those 
leading effects for several mesonic two-point functions (pseudoscalar, scalar singlet, axial 
and vector). The interesting feature of this regime is that only the O(a^) corrections in the 
chiral Lagrangian contribute. There is no proliferation of unknown couplings: apart from 
the continuum leading order couplings S and F, only an extra constant C2 appears. 

We have computed the leading O(a^) corrections to the PCAC quark mass and expressed 
the correlators as a function of the dimensionless variables Zm = rnpcAc^V and 2^ = fi'EV. 
We have adopted the so-called twisted basis, where isospin breaking for non-zero twist angle 
explicitly shows up. 

The final formulae we quote are valid for an arbitrary twist angle, and hence reproduce 
also the untwisted Wilson case considered in [381 ES]! for a; = (equivalent to /i = 0). A 
particularly interesting setup is maximal twist u = tt/2, here defined by mpcAc = 0, where 
automatic 0(a) improvement occurs ^U\. The numerical investigations performed in section 



IV D| suggest that, like for untwisted Wilson fermions, for typical lattice parameters adopted 
in present-day simulations, the O(a^) corrections remain at the few percent level. This result 
supports the possibility to extract low-energy couplings with twisted mass Wilson fermions 
simulations in the e-regime with controlled systematic errors. Notice, however, that it is 
not possible to predict a priori in which particular regime (GSM, Aoki) one actually has 
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FIG. 3: The axial vector correlator with a = 1, 2 (normalized by —F'^/T) as a function of z 



performed a simulation, and a scaling study is advocated. 

We finally remark that a determination of the twist angle in the e-regime might be 
numerically not so easy. To determine the actual value of u with a reasonable accuracy the 
statistical and systematic uncertainties on the PCAC mass have to be ideally much smaller 
than the value of the twisted mass /i. In the e-regime, where the value of the twisted mass 
could become comparable to the errors associated with the PCAC mass, this task could 
become extremely difficult. This fact could induce large uncertainties in the determination 
of the twist angle. A realistic procedure is to determine the twist angle from p-regime 
simulations [5Sl EI] where it is possible to keep the uncertainty on the twist angle well 
under control. The bare parameters tuned in the p-regime are then used in the e regime. 
We expect that this choice will induce 0(a) corrections to the PCAC mass, which will not 
spoil automatic 0(a) improvement and the validity of the formulae given in this paper. 
An alternative procedure is to use correlation functions which are u independent at the 



classical level. As we have seen in sect. IV C (cf. eqs. (54) and (62)) these particular linear 



combinations are u independent also at NLO in the GSM* regime and additionally they are 
free from O(a^) corrections. This is an alternative procedure to analyze lattice data which 
is free from the uncertainties stemming from the determination of the twist angle and free 
from the O(a^) corrections affecting the standard correlation functions. 
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Appendix A: Selected formulae and definitions for the epsilon regime 



In this Appendix we briefly summarize formulae which are relevant for the computation 
of correlation functions in the e-regime of chiral perturbation theory. For more details 
the reader can refer to [381 EH]- The pseudo Nambu-Goldstone bosons propagators for the 
nonzero modes is defined as 



ipx 



(Al) 



Following Refs. 



we define 



G(0) 



/9i 



V 

00 



V 



(A2) 
(A3) 



where /3i and /cqo are finite dimensionless shape coefficients which depend on the geometry 
and can be evaluated numerically. 

The integral over spatial components of the propagators G{x) yields the parabolic function 

2/ 12 



/ 



d^xG{x) = Thi 





T 




t 


(^) 


~ 2 


[( 


T 



(A4) 



Another quantity which appears frequently is the quark condensate at one loop, which 
for Nf = 2 is given by [S2] 

3 /3i 



^cff 



S I 



2F2 ^) ■ 



(A5) 



In the calculation of the NLO continuum correlators the measure factor in eq. (25) is 
needed at ©(e^) [H [52]: 



B{0 



4 1 
3F2y2 



c/^a;Tr(^'^T"^^T*^ 



(A6) 
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Appendix B: Continuum correlators in the twisted basis 



In this appendix we summarize continuum formula for two-point functions in the twisted 
basis. For simphcity we introduce the function 



L(x) 



xh{2x) 

The results are given in terms of the variables 

Zf, = /iSV, Zm = mpcAcSV, 



(Bl) 



(B2) 



The pure Wilson case corresponds to z^ = 0, z = Zm, while the maximal twist is verified 
for Zm = 0, z = Zfj,. Each correlator consists of a constant part and on a time-dependent 



contribution; the time-dependence is represented by the function hi defined in eq. (A4). 



1. Scalar and pseudoscalar correlators 

We write the continuum scalar singlet correlator as 

C's%At)=L'Cl'^ + a'^s'Thit/T) 

while for the pseudoscalar correlator we have 

CfpJt) = L^C^J' + afTh{t/T) 

The coefficients C50, aso,Cf^ and ap are given by 



a 



00 



C 



s 
11,22 



F2 



Z^ 
-v2 



1 



+ 



Liz) 



1 



2(^^ - zl 



2 ^(^eff) 



a 



cf = 


y2 


[z^ 

_Z^ 


11,22 




1 - 


P ~ 


F2 


^33 




m 


ttp = 


F2 





Liz. 



cS) 



1 + 



L{z) 



+ 



2 

Liz 



1 + 



2iz^. - zl 



(B3) 
(B4) 

(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 



The subscript "eff" in Zes indicates that S must be substituted by the one-loop corrected 



quark condensate given in eq. (A5). 
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2. Axial and vector correlators 



The continuum vector and axial correlators are given by 



^AAyV,ct^ 



X 



y) 



a' 



ab 

Ay 



T 



+ ^Ptvkoo-^ltvhi{t/T). 



(Bll) 



Here Pi and fcoo are familiar shape factors which depend on the geometry of the space-time 
volume. They are defined in eqs. (A2),(A3). 
Explicit results for the coefficients are: 

9 ^2 _ ^2 



a, 



11,22 



+ 



2/3i 



2 



+ 



Liz] 



a 



33 



2/3i 



(1 - L(;.eff)) + ^ (1 - L{z)) 



oll,22 

Pa — 




PT = 


2L{z) 


11,22 

7a = 




7f = 


2z^L{ 


11,22 

ay' 


= F^ 



+ 



-L{z) 



+ 



z^ - 

-^^Liz, 

Z'^ 



cflf. 



+ 



+ 



+ 



^2 -^2 



a 



33 
V 



F^L(z,s) + ^L(z) 



11,22 
V 



-Liz] 



11,22 

7y 
Iv 



2(1 
0. 



Liz)) 



(B12) 



(B13) 

(B14) 

(B15) 
(B16) 
(B17) 

(B18) 



(B19) 

(B20) 

(B21) 
(B22) 
(B23) 



Also in this case, the subscript "eff" in Zes indicates that S must be substituted by the 
one- loop corrected quark condensate Hes given in eq. (A5). 



Appendix C: Notations for comparison with ref. 

In ref. [38] a slightly different notation has been used in comparison with this work. Here 
we give a short summary of the main differences. In ref. |3S] the variables z and Z2 have 



been introduced. They correspond respectively to Zm/2 (cf. eq. (32)) and — p (cf. eq. (30)) 
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In the following for the z variables we will use the same notation adopted in the main 
text given in eqs. (46) and (47). The Bessel functions /n(x) translate to the function 
used in [SB] as 

(^/2)" 



Ar(n + i/2)^"- 

As examples we give here few ratios which typically appear in correlation functions 

/i ~ 3 Xi ~ Xi ' Ji ~ 15 Xi ' 
To conclude we find useful we write the following correlation function 



(CI) 



(C2) 



with the notations of ref. 



We have 
X2 



Cpp,ct(^) + 



11 

PP,a? 



^cff 



Xi ^ F2 



1 - 



1X2 
3X1 



T 

Z3 



hi{t/T) 



(C3) 



(C4) 



and 



^11 




( 32 



+ 



3^2 



X2 



(C5) 



where X„ = Xn{2z) with z defined in eq. (47). 



Appendix D: Example of group integrals with isospin breaking 



In this appendix we show one of the procedures we have used to perform the zero- 
modes integrals with isospin breaking integrands. We use the example of calculating the 



expectation value {dSa'^), which is part of the correction (28) we are interested in. To 
simplify the calculation we restrict ourselves to the case of maximal twist with m = 0, 
which nevertheless demonstrates the main idea. In this case the integral we have to perform 



reads (c.f. eqs. (29) and (31)) 



(SS^^) = ^ / d[f/o](Tr(f/o + t/t))V^^^^[-^(^-^o^)l. 
loZo J 



(Dl) 



Performing the change of variables Uo — )■ ia^Uo (which is equivalent to the rotation (|5| with 
uo = 71/2) this translates into 



(SSa^) = -T^ I d[t/o](Tr(a3[t/o-f/t]))^e^^^[^o+^or 



16Zn 



If we use the parametrization 



Uq = exp[i(pn ■ a/2] 



(D2) 



(D3) 
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(0 G [0, 27r] , \n\ = 1) for the constant mode, the Haar measure reads 

d[Uo] = ^d^dfisin^^, 

with dQ = dQ{n) being the measure of the 2-sphere, and the integral turns into 

1 



/ d[f/o]e' 



2zu cos f 



sm —Tin 



(D4) 



(D5) 



The factor n| in the integrand is a remnant of isospin breaking. Without isospin breaking 
the integrals one faces involve functions of (p only, at least if various completeness relations 
for the group generators are used [S2j- In this case the integration J dQ gives a trivial 
factor 4-7r and the remaining integral over leads to expressions involving modified Bessel 
functions. In order to make contact to these known integrals we write 



dQnl = - 

and the integral we are interested in turns into 



dQ, 



Pi 
3Zn 



d[Uo]e 



2Za COS 



(D6) 



(D7) 



For the remaining integration over we write sin^(0/2) = 1 — cos^(0/2) = 1 — (Trf/o)^/4, 
and we finally obtain 



P 

3 



(l)phys - ^((Trf/o)^)phys 



(D8) 



This is our desired result: The right hand side involves familiar integrals with the standard 
Boltzmann weight, c.f. ([33|) (as indicated by the subscript "phys"). A useful collection of 



relevant integrals is given in appendix B of Ref. [39], which can be used to express (D8) as 

1 /2(2^^) 



(D9) 



The same steps can be carried out in the calculation of correlators {0^0^ . The integrand 
will be a product f{(j))p{n3), where / is a function of only and p denotes a polynomial in 



n^. Generalizing eq. (D6), the integration over S is trivial and gives a simple factor c„. 



dQp{n3) 
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dQ. 



(DIO) 
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The remaining integral can be expressed in terms of familiar integrals known from continuum 
e-regime calculations without isospin breaking. 
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